Green's Functions from Quantum Cluster Algorithms [] 

'66 : 

■ R. Brower a,c , S. ChandrasekharanPf 3 ' and U.-J. Wiese c 



-a 



a Department of Physics 



! Boston University 

Boston MA 02215, USA 

t— i ■ 

a T-8, Theoretical Division 
' c| Los Alamos National Laboratoy 

g i Los Alamos, NM 87545 

a: 

c Center for Theoretical Physics 
Laboratory for Nuclear Science and Department of Physics 
^ ! Massachusetts Institute of Technology 

i Cambridge MA 02139, USA 

a 

_L ■ Preprints 

MIT-CTP 2686 
LA-UR 98-91 

^ ; PACS: 02.70.Lq, 05.30.-d, 05.50. +q, 11.15.Ha 

Keywords: Quantum Spins, Clusters, Green's Functions, Improved Estimators 

cn 

O 
O 

o 

00 

ON . Abstract 

We show that cluster algorithms for quantum models have a mean- 
ing independent of the basis chosen to construct them. Using this idea, 
we propose a new method for measuring with little effort a whole class 
of Green's functions, once a cluster algorithm for the partition function 
has been constructed. To explain the idea, we consider the quantum XY 
model and compute its two point Green's function in various ways, show- 
^ ■ ing that all of them are equivalent. We also provide numerical evidence 

confirming the analytic arguments. Similar techniques are applicable to 
other models. In particular, in the recently constructed quantum link 
models, the new technique allows us to construct improved estimators 
for Wilson loops and may lead to a very precise determination of the 
glueball spectrum. 
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1. Introduction 



Cluster Algorithms have had a tremendous success in numerically solving some 
interesting field theoretic models 0]. Such algorithms have the potential to dras- 
tically reduce the problem of critical slowing down, and perhaps even eliminate it 
completely. The Swendsen-Wang || algorithm was the first cluster algorithm to 
be invented. It has been extended to both classical Q and quantum spin systems 
P, 0] with continuous symmetries. On the other hand, the success of cluster al- 
gorithms has been limited to spin models without frustration. Further, it has not 
been possible to construct an efficient cluster algorithm for Wilson's formulation of 
lattice gauge theory. However, a new formulation of gauge theories, proposed re- 
cently, provides a possible alternative. The corresponding quantum link models |7j 
are a quantum version of Wilson's gauge theory, in the same sense as the quantum 
Heisenberg model is a quantum version of the 2-d 0(3) model. A cluster algorithm 
for the simplest U(l) quantum link model has recently been constructed While 
in quantum spin models the clusters are closed loops, in quantum link models they 
are world-sheets of electric flux strings. It is expected that similar algorithms exist 
for other gauge groups. Thus, it appears that cluster algorithms can also play an 
important role in the numerical solution of gauge theories. 

Updating a system with an efficient cluster algorithm accelerates its evolution 
through configuration space. However, in quantum models this does not automati- 
cally imply that one can efficiently measure various Green's functions. For example, 
one is interested in creation and annihilation events described by operators that are 
not diagonal in the basis in which the path integral is constructed. The measurement 
of such "quantum flip" operators appears to be very difficult. In fact one appears 
to be restricted to a very limited class of Green's functions built from the diagonal 
operators. For example, let us consider the quantum XY model expressed in terms 
of the Pauli matrices a l x at the lattice sites x. When we work in the basis where a\ 
is diagonal, Green's functions for products of a x operators can be measured easily. 
On the other hand, it appears that Green's functions involving a x or a x are difficult 
- if not impossible — to measure. A measurement of a\ or a x in the path inte- 
gral described by eigenstates of a x necessarily "flips" the state of the corresponding 
quantum spin. In other words, Green's functions of a x or a x receive contributions 
from superselection sectors of the theory, which are not sampled by a conventional 
Monte Carlo algorithm which simply simulates the partition function in a particular 
basis. 

For quantum cluster algorithms, there is a remarkable way out of this impasse. 
By virtue of the auxiliary percolation process used in a quantum cluster Monte 
Carlo algorithm, the requisite information for all quantum Green's functions can in 
fact be extracted. In this paper, we formulate a method, which we call "quantum 
flip" measurements, for computing non-trivial Green's functions from a quantum 
spin cluster simulation. Our construction of the "quantum flip" measurement pro- 
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cedure is based on a generalization of the classical Fortuin-Kasteleyn || mapping to 
quantum clusters and on the observation that a basis-independent description exists 
for the resulting quantum random cluster model. For the Swendsen-Wang clus- 
ter algorithm, it is well known that the Fortuin-Kasteleyn mapping demonstrates 
the equivalence between the partition function of the Ising or Potts model with 
the partition function of an associated random cluster model. This step, which is 
sometimes ignored in discussions of cluster algorithms, leads to a purely geometric 
procedure to measure spin correlations. For quantum spin models our generalized 
Fortuin-Kasteleyn map naturally also leads to a systematic and rigorous method to 
measure any operator or Green's function. The quantum Green's functions involv- 
ing operators non-diagonal in a particular cluster update operation are related to 
cuts (or "flips") which cleave the original cluster into disjoint sub-clusters. Based 
on this picture, it is possible to apply similar ideas to more complicated systems 
such as gauge theories expressed as quantum link models. In qualitative terms a 
similar suggestion was made in ref.[|10j. Our "quantum flip" algorithm gives a sim- 
ple explanation of the idea and provides a general approach to measure any Green's 
function. 

The article is organized as follows. In section 2, we reformulate the quantum 
spin model as a model that describes the dynamics of clusters with internal quantum 
degrees of freedom. We use the quantum XY model as an example. The reformula- 
tion demonstrates the basis-independence of the quantum clusters. In section 3 we 
show how one can measure any operator using the cluster algorithm. In section 4, 
we illustrate how the observations of sections 2 and 3 are realized in the quantum 
XY model, by constructing two cluster algorithms which are related by a change of 
basis. We show why the cluster properties are independent of the basis, and provide 
numerical evidence confirming the analytic results. In section 5 we suggest how to 
use similar ideas in the U(l) quantum link model to obtain improved estimators for 
Wilson loops and present our conclusions. 

2. Basis-Independence of Clusters 

In this section, we show that constructing a cluster algorithm for a quantum spin 
model is equivalent to rewriting the path integral in terms of geometrical objects, 
which we call quantum clusters, with internal quantum degrees of freedom, in this 
case quantum spins. We then argue that the dynamics within each cluster can be 
written in a form which is invariant under basis transformation of the quantum 
spin Hilbert space. In the next section, we discuss how these observations lead to a 
method for measuring a whole class of operators. In order to illustrate the idea we 
consider the two-dimensional quantum XY model. Most of the observations extend 
trivially to other models including quantum link models. 
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The quantum XY model is described by the Hamilton operator 

H = J2 h ^ ; h *,fi = g { a x ff l+pi + a l a l+t) ■ (!) 

Here ft = l,2,...,d, denotes spatial directions and x denotes the position on an L d 
lattice. The partition function is given by 

N 

Z = Trexp(-/3#) = Tr [] exp(-e#). (2) 

i=l 

In the last step we have introduced the Suzuki- Trotter discretization such that eN = 
p. It should be noted that it is not necessary to discretize time. For quantum 
systems in a discrete basis the path integral is well-defined in continuous time, and 



can be simulated directly in the time-continuum ||11|| . In this paper we work at 
discrete time, but all results have a direct continuous-time version. In the above 
form the partition function is not yet accessible by numerical simulations. In fact, 
the term exp(— eH) needs to be expanded further. The basic idea is to write the 
Hamiltonian as a sum of interactions arranged in a checker board pattern such that 
a minimal number of spins interact in a single time step. For example, in d = 2, one 
decomposes 

H = H 1 + H 2 + H 3 + # 4 . (3) 
This divides the Hamiltonian into four terms in a checker board pattern with 

x&(2m.n) x(^(m,2n) i6(2m+l,n) xG(m,2n+l) 

(4) 

Every Hi contains a sum of commuting operators, each of which represents a two-spin 
interaction. Thus exp(— eif;) can be computed easily, and evolves all the spin states 
through a single time step. The partition function is then defined on a hypercubic 
lattice with AN time slices. It is equal to the trace over the Hilbert space of a product 
of 2NL 2 transfer matrices connecting pairs of spins. This is because the spins 
interact in pairs in a checker board pattern depending on which of the exp(— eHj) 
is active in a particular time step. There are V/2 two-spin interactions, where 
V = 4NL 2 is the volume of the lattice. In the d — 1 case, one only needs two 
terms in the expansion of the Hamiltonian, i.e., H — H\ + Hi. Then there are 2iV 
time slices and the partition function is a trace of a product constructed from NL 
transfer matrices, connecting pairs of spins. A typical lattice for d — 1 is shown in 
figure 1. Generalization to higher dimensions is straight forward. 

The transfer matrix that evolves a pair of spins through a single time step is 
denoted by an operator T. Since T represents a two-spin interaction, it is equivalent 
to a 4 x 4 plaquette matrix. Denoting the basis for the spins by \s), the operator T 
can be expanded as 

T= l s 3S 4 )'?"(si, s 2 , s 3 , s 4 )(sis 2 | (5) 

Sl,S2,S3,S4 
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where s±, s 2 are the spins in one time slice and S3, s 4 are the same spins in the next 
time slice. For example, if we choose the basis to be eigenstates of a 3 ,, the non-zero 
matrix elements of T are given by 

n+ + ++) =r( ) =1, (6) 

T(+-+-) =T(- + -+) =cosh(eJ), 
T(+--+) =T(- + +-) =sinh(eJ). 

Here ± represent the eigenstates of a 3 with eigenvalues ±1. This leads to the path 
integral form for the partition function, 

V/2 

Z = Tr(YlT p )=Y l W[si\, (7) 

p=l <S 

where of course the product over operators T p for each of the V/2 plaquettes must be 
appropriately "time" ordered. In the last step, we have replaced the trace by a sum 
over spin configurations S, which is a collection of spins Sj located at the space-time 
points of a (d + l)-dimensional lattice, weighted by the weight W. The associated 
probability distribution of the spin configuration is given by V[si] = W[si]/Z. The 
weight function is a product over the four-spin plaquette interactions, which takes 
the form 

W[s t ] = T(S 1 , S 2 , S3, S 4 )T(S 4 , S 5 , S 6 , S 7 )...T(sy_ 3 , Sy_ 2 , 8 V -\, S V ). (8) 

Due to the trace structure of the partition function, each spin state Sj occurs in 
two different T operators. The sum over a particular spin Sj makes the partition 
function independent of the basis that is used to describe that spin state. We will see 
that the basis-independence of the partition function extends to the clusters. One 
can rewrite it completely in terms of clusters without referring to the basis of the 
Hilbert space. This observation also leads to similar definitions for the expectation 
values of operators. 

We now proceed with the construction of a cluster algorithm and the Fortuin- 
Kasteleyn map onto its associated "random cluster" model. The construction of a 
cluster algorithm is equivalent to rewriting T(si, s 2 , s 3 , s 4 ) as a sum over products 
of simpler tensors For example, we can expand 

T(si, s 2 , s 3 , s 4 ) = ao-MoOi, s 2, s 3 , s 4 ) + aiMi(si, s 2 , s 3 , s 4 ) + a 2 M 2 (s 1 , s 2 , s 3 , s 4 ), 

(9) 

as a sum of three terms, where we selected a minimal tensor basis, 

M = 5 S1 , S3 5 S 2,S4, Mi = 5 SuSi 5 S2 ^, M 2 = o"s llS2 o-s3, S 4' ( 10 ) 
with the respective weights uniquely determined to be 

a = |(1 + exp(-eJ)), a x = |(1 - exp(-eJ)), a 2 = i(exp(eJ) - 1). (11) 
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Thus each tensor T in the partition function is equivalent to a sum of three tensors, 
M. n , labeled by a plaquette percolation state, n p = 0, 1, 2. This allows us to rewrite 
the partition function as a sum over 3 V ^ 2 terms since there are V/2 Ts in the partition 
function. By choosing one of the three terms for each plaquette, one is led to a 
joint probability distribution, V[si, n p ], where n p labels the percolation state of each 
plaquette. Equivalently, since we have chosen a tensor decomposition in which each 
term is given by a matrix connecting two adjacent spins, each of the 3 y / 2 terms can be 
represented by a graph G with connected components (i.e. clusters) Cj, i — 1, Nq. 
A simple example of a graph in a 1-d model model is shown in figure 1. 




Figure 1: A typical graph in 1-d with L = A, and N = 2: Since H = Hi + H 2 , 
one has 2N time-slices and V = 2NL. Every white square represents a two-spin 
interaction. The graph G shown here splits the lattice into three clusters, represented 
by a solid line, a broken line and a dotted line. It is easy to read off w{G) = a\a\a\. 



The Fortuin-Kasteleyn map is based on the simple observation that the joint 
distribution , V[si;n p ], has two marginal distributions: One, obtained by summing 
over the percolation states, is the original spin probability distribution, 

V[si] = £Ph;rg, (12) 

{n p } 

while the other marginal distribution, obtained by summing over the spin variables, 

V[G]= J2V[ Sl ; n p ], (13) 

5 
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is a new "random cluster model" . The notation on the left-hand side acknowledges 
the fact that in our application each distribution of integer percolation states n p 
uniquely corresponds to a graph G. 

It should also be noted that this mapping can be cast in purely operator terms 
for the quantum model. Going back to the operator form of the partition func- 
tion, we have a plaquette operator, T, in eq.@ for each four-spin matrix element, 
T(si, S2, S3, S4). Obviously, the operator T can be expressed as a sum of simpler op- 
erators, J2n a nM n with matrix elements that reproduce the above expansion eq.@ 
for T(si, s 2 , s 3 , s 4 ). The contribution to the partition function for a particular graph 
G takes the form of a trace Tr[M(j], where the operator Mq is the appropriate time- 
ordered product of factors M np and the full partition function can be rewritten as 

Z = J2w(G)Tr[M G ]. (14) 

G 

Thus we have defined the probability distribution, 

V[G] = i^(G)Tr[M G ], (15) 

for the associated quantum random cluster model. 

The cluster algorithm is a way to realize these dynamics in a Monte Carlo 
method. In a particular basis it reproduces a graph G with the clusters Ci(i = 
1,...,Ng) with the weight w(G) and internal spin orientations Tr[M<?]. However, 
since w(G) and Tr[Mc] are basis-independent, it is clear that algorithms can be 
constructed in more than one basis. Such algorithms produce clusters with exactly 
the same geometrical properties. We will illustrate this in the case of the quantum 
XY model in section 4. It is possible that the efficiency of algorithms working in 
different bases are different since the cluster growth dynamics would be different. 
Moreover, it is important to realize that in a new basis one may choose new cluster 



decompositions defined by eq.(|y) |L2| which are not equivalent to a basis change. 
In this case, the cluster properties of two random cluster models are in general 
inequivalent. 

The above picture is a generalization of the well-known Fortuin-Kasteleyn map- 
ping H to quantum models. Such a quantum mapping was also considered in 



another more technical paper ||13fl . Here we present a new way to understand the 
mapping. For example, we have shown that the quantum XY model is equivalent 
to a model that involves the dynamics of geometrical objects (in the present case 
closed loops), with internal quantum degrees of freedom (in this case spins). Sum- 
ming over the spin degrees of freedom defines a partition function that describes 
the dynamics of the clusters independent of the basis used for the spin states. This 
is at the heart of the Fortuin-Kasteleyn mapping for the Potts model. As we will 
show in the next section, the basis-independence of the clusters brings into focus 
some interesting ideas on how to measure general Green's functions. In particular, 
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it shows why the idea of ref.|]IU| works and provides a systematic way to measure 
any operator. 



3. Measurement of Green's Functions 



The basis-independence of the clusters suggests that it must be possible to mea- 
sure any operator with a given cluster algorithm. This is indeed true. For example, 
in the XY model it is possible to measure the <j].<Jy correlators even when one uses 
the cluster algorithm in the basis where a 3 is diagonal. The fact that such mea- 
surements are possible has recently been suggested in ref.|H|]. Based on the above 



ideas, we explain how this is possible. In fact, our approach extends to all oper- 
ators, although the measurement of the above correlator seems to be the simplest 
illustration of our approach. 

The expectation value of an operator O x , acting at the site x, is defined by 

(O x ) = ±Tr[O x exp(-pH)]. (16) 
Rewriting the above expression one obtains 

(O x ) = ^ w gTt[O x M g ]. (17) 

L G 

Let us consider the steps that lead to the numerator of eq.(|17]). The numerator of 
eq.flTBT) can be written as a product of transfer matrices T as described in the previous 
section. However, now there is an additional operator O x acting on a specific site in 
a particular time-slice t. Then there are two transfer matrices adjacent to this site. 
The matrix T~ connects the site with the previous time-slice, and T + connects it 
with the next time-slice. The operator O x is sandwiched between the two. In the 
next step each T is expanded into three terms and the whole expression is again 
written as a sum of 3 V ^ 2 terms. Each of the 3 V ^ 2 terms represents a graph G with 
clusters C iy (i = 1,...,N G ). The weights w{G) are the same as before. However, 
due to the presence of the operator O x the previous expression is modified. It now 
contains the matrix O x appropriately sandwiched between the tensors A4q, -Mi, or 
M.2 coming from T + and T~ . In the example of the quantum XY model, where the 
clusters are closed loops and M G is a product of 5 ss i or (cr 1 ).^, one simply inserts 
the matrix (O x ) ss i at the right place. When the cluster turns back in time, the 
appropriate matrix to be inserted is (0^) s>s >. Such subtleties must be kept in mind. 

It is easy to find the right prescription to measure any operator using the cluster 
algorithm. The cluster algorithm produces the graph G with the weight, 

V[G] = i w (G)Tr[M G ]. (18) 
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Thus we get the formula 

_ Tr[O x M G ] 
{ ° x) ~ ( Tr [Mq] } - (19) 

To evaluate this expression one must perform all possible cluster flips, which can 
be done very efficiently using an improved estimator. As we will see in the next 
section, the above formula can be easily extended to measure any two-point Green's 
function in the quantum XY model. 



4. Two-Point Green's functions 

Here we will demonstrate the power of the above ideas by measuring various two- 
point Green's functions in the quantum XY model, using the cluster algorithm in the 
a 3 basis. We show how the dynamics of the clusters allows us to prove analytically 
that the 0(2) symmetry is maintained. We also construct a cluster algorithm in a 
different basis, and show how the same features emerge in the new basis. 

Let us consider a cluster algorithm obtained by decomposing the transfer matrix 
as given in eq.(|9|). The rules for cluster growth are given in figure 2. One can 

s • • s 



l l 

s • o -s 



P3 1-P3 



-s s 
s #^.~>o -s 

Figure 2: The rules for creating the bonds that define the graph G: The sites are 
described by filled circles and the time axis is vertical. The spin states in the a 3 basis 
are denoted by s, which can be ±. The probabilities for choosing the bonds are pi = 
[l-exp(-eJ)]/2, p 2 = [exp(eJ)-l]/[2cosh(eJ)], andp 3 = [exp(eJ)-l]/[2 sinh(eJ)]. 

easily derive some properties of the clusters. They are closed loops, since each 
site is connected to two bonds, one for an interaction forward, and the other for 
an interaction backward in time. Further, the bonds connect spins of the same 
orientation when the cluster grows forward or backward in time. However, when 
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they connect spins in the same time slice, they connect opposite spins. Globally, this 
leads to an interesting property of the clusters. If the clusters have positive spins 
when they grow forward in time, they have negative spins when they turn back and 
grow backwards in time. This property of the clusters is important for proving that 
the correlators (&l<Jy) and (cr^cx 2 ) are equal. 

Let us measure the correlator (a^cr*) with the above cluster algorithm. Note that 
since a\ changes the spin state, it is a priori unclear how to measure {p\a^) , given 
just the spin configuration. In fact, a spin configuration that contributes to the 
correlator will never be generated in the cluster algorithm. However, eq.(|T9|) tells us 
to measure Tr[al<TyMG}/Tr[Mc} for each graph generated by the cluster algorithm.^ 
Each graph G represents a set of clusters Cj (i = 1, N G ). Since M G is a product 
of 8 ss i or a] s ,, the trace contracts the spins on each cluster separately. Thus Tr[Mg] 
becomes a product of traces for each cluster, i.e., 

N G 

Tr[M G ] = nTr[M c J, (20) 

where N G is the number of clusters in G. Further, given one spin in the cluster, 
the orientations of all other spins in the cluster are determined. This means that 
Tr[McJ =2. In other words, TrfM^] measures the number of configurations ob- 
tained from all possible flips of the clusters. Since each cluster contributes two 
configurations, Tr[Mc] = 2 N ° . On the other hand, TrfcrJcr^MG-] is the same, except 
that at the sites x and y the spins are now flipped. Clearly, if x and y were in dif- 
ferent clusters, this flip would be inconsistent going around the closed loop, which 
means that Trfa^McJ = 0. However, when x and y are both in the same cluster, the 
contribution is exactly the same as Tr[M G ]. Thus, each time a cluster is produced 
that includes both x and y, we get a contribution 1 to the measurement of (cr^cr*). 
Due to the translational invariance of the Hamiltonian, it is possible to consider 
any pair of points x, y within a given cluster and count it as a contribution to the 
correlator. In this way, it is possible to obtain the complete two-point correlation 
function (o^cr*). 

What about the (cr 2 cr 2 ) correlator? Since the model is 0(2) invariant, it must 
be equal to (o^a*). This is indeed the case, although the proof is a bit subtle. 
While an insertion of a 1 flips the spin, an insertion of a 2 flips the spin along with 
a multiplication by a factor ±z. This means that, when x and y are in different 
clusters, Tr[cr 2 M c .] is zero as before. However, if they are in the same cluster, one 
must figure out what the factors ±i contribute. Here we invoke a special property of 
the clusters to show that there is always exactly one i and one —i factor coming from 
the two a 2 insertions. The point is that if x and y are in the branch of the cluster 
that grows in one particular time direction, the spin flips from +(— ) to — (+) at the 
first cr 2 insertion and from — (+) to +(— ) at the second insertion. This gives one 

3 In the interest of simplicity we ignore the fact that when the operators a\ and a y are on 
different "time" slice the product a^alMc should, of course, be appropriately "time" ordered. 
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factor i and one —%. On the other hand, if x is in a branch growing in a particular 
time direction and y is in a branch growing in the opposite direction, the spin flips 
from +(— ) to — (+) at the first and from +(— ) to — (+) at the second a 2 insertion. 
However, since the cluster is growing backward in time at one of the insertions, the 
actual matrix element is that of (c 2 ) T which again gives one i and one —i factor. 
Thus, if x and y are in the same cluster, the contributions to (cr^a*) and (cr^cr 2 ) are 
the same from each graph generated by the cluster algorithm. 

It is straightforward to measure (o"^cr~) or (o~~o-y) in the same way, and one 
gets half of the contribution of (cx^cr^) for each. Thus, it is possible to measure 
any two-point Green's function directly using formula eq. (|19D . In fact, one can 
easily use it also to measure higher n-point Green's functions, which, however, get 
contributions from more than one cluster. It is interesting that there is a cluster- 
notion of "connected" and "disconnected" pieces to the four-point Green's function. 
The connected piece comes entirely from within a cluster, while the disconnected 
piece comes from two separate clusters. 

We now consider a change of basis and show the basis-independence of the clus- 
ters using a concrete example. One can consider a basis change given by 

k> = E( e_ia * 7r/4 )— t K>- (2i) 

4 

This basis change can be applied to each of the spin states. In terms of the new 
basis, the transfer matrix given in eq.(|5]) changes to 

Sl,S 2 ,S 3 ,S4 s'^s'^s'^ 

144) T(si,s 2 ,s 3 ,s A ) (44l(e wl7r/4 ) s ' lSl (e^ l7r / 4 ) 4s2 (22) 

The transformed tensor can be written as 

T'isx, s 2 , s 3 , s 4 ) = a A^' (si, s 2 , s 3 , s 4 ) + aiM[(si, s 2 , s 3 , s 4 ) + a 2 M' 2 (s 1 , s 2 , s 3 , s 4 ) 

(23) 

where 

= ^Sl,S3^S2,S4) -M-l = ^Si,S4^S2,S3) ^M- 2 = ^ Si ,S2 ^ S3 ,S4 ■ (^4) 

Notice that the change of basis has only effected A4 2 , leaving A4 and A4i un- 
changed. It is clear that all the arguments presented in the previous section are still 
applicable here. In fact, one can construct a cluster algorithm for the quantum XY 
model in this basis, which is the same as an algorithm for a quantum XZ model in 
the <t 3 basis. The cluster rules are given in figure 3. 

Also in the new basis the clusters are closed loops. However, now all the spins 
in a given cluster have the same orientation. Since now we are dealing with the XZ 
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Figure 3: The rules for creating the bonds that define a graph G in the XY model in 
the basis in which a 2 is diagonal: The sites are described by filled circles and the time 
axis is vertical. The spin states in the a 2 basis are denoted by s, which takes values ±. 
The probabilities for choosing the bonds are given byp\ = [1— exp(— eJ)]/[exp(eJ)+l] 
and p 2 = [1 + exp(— eJ)]/[exp(eJ) + 1]. 

model, the (cr^cr^) correlator must be the same as the {or 2 o- 2 ) correlator in the XY 
model. It is easy to see that, if x and y belong to different clusters, Trfer^MeJ is 
zero. However, when they belong to the same cluster, since all spins within a cluster 
have the same orientation, Tr [cr^o^MeJ = Tr[McJ. Thus, one gets a contribution 1 
to (cr^er^). This is exactly what we had in the previous basis for (cr 2 a 2 ). Note that 
in the previous basis one cannot show that (cr^o^) gets a contribution 1 whenever 
x and y belong to the same cluster, since opposite spins can be present in the same 
cluster. 

We have obtained numerical results that confirm the above conclusions. We have 
updated a 10 x 10 lattice with N = 100 at (3 J = 1.0 using both the XY and the XZ 
algorithm. The above arguments imply that both algorithms must yield identical 
cluster properties. The cluster size distribution has been measured with a single 
cluster algorithm. Figure 4 shows the histograms for the cluster size distributions 
obtained with the two algorithms. We see that they are identical within statistical 
errors. In table 1 we give the numerical values for the first few moments of the 
cluster size distributions of figure 4. 
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Figure 4: The histograms of cluster size distribution for the two algorithms discussed 
in the text: Both algorithms were run with single cluster updates. The data sets 
contain 10000 clusters each. The results of the XY algorithm are given by the shaded 
and that of the XZ algorithm are given by the unshaded histogram. 

5. Conclusions 

In this article we have developed a new method to evaluate any Green's function 
in quantum spin models using cluster algorithms. The method is based on the fact 
that clusters have a meaning independent of the basis. In other words, the construc- 
tion of a cluster algorithm implies that one can rewrite the quantum spin model as a 
random cluster model in a basis-independent way. This is well-known in the case of 
the Potts model using the Fortuin-Kasteleyn mapping |9j . Here we have presented a 
simple argument showing that such a mapping also exists for quantum spin models. 
This mapping provides a new method to measure any operator once a cluster algo- 
rithm has been constructed. In this method a measurement of a certain operators 
reduces to the computation of traces in the Hilbert space for each individual cluster. 
In the quantum XY model, this turns out to be quite simple. The method applies 
to other models as well, and is useful for measuring quantities which could not be 
measured using standard methods. 

The insight gained here extends easily to quantum link models, which were re- 
ferred to in the introduction. Since these models are gauge invariant, the operators 
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Algorithm 


(C) 






XY 


9.750(120) x 10 3 


1.790(33) x 10 8 


3.941(96) x 10 i2 


XZ 


9.818(119) x 10 3 


1.807(34) x 10 s 


3.990(95) x 10 12 



Table 1: The first few moments of the cluster size for the two algorithms: The 
simulations were performed on a 10 x 10 lattice at (3J = 1.0 and N = 100. This 
clearly demonstrates the basis-independence of the clusters. 



of interest are closed Wilson loops. It seems difficult to construct cluster algorithms 
in the basis where the Wilson loops can be measured easily. For example, in the 
£7(1) quantum link model, a cluster algorithm has been constructed in an electric 
flux basis. The clusters resemble world-sheets of electric flux strings. When Gauss' 
law is imposed, the flux strings form closed loops and hence sweep out closed world- 
sheets. Although the Wilson loop operators are not diagonal in the electric flux 
basis, they can be computed using the method proposed in this article. The com- 
putation requires understanding the topology of the closed world-sheets generated 
by the Monte Carlo algorithm. For example, only if the cluster surface separates 
into two pieces when it is cut along the links belonging to the Wilson loop, there is 
a non-zero contribution to the operator. It is easy to show that Wilson loops get 
contributions from only one cluster at a time. At present, this is used to investi- 
gate the physics of the U(l) quantum link model. The results of this study will be 
published elsewhere 0. 

The connection with the random cluster model should be understood in more 
detail in the context of quantum cluster algorithms. One intriguing possibility is 
a solution of the negative sign problem, at least in a class of models. Perhaps the 
random cluster model suggests a natural basis to avoid negative signs. The new 
insight may also serve as a guide for constructing cluster algorithms for quantum 
link models with non-Abelian gauge groups. 
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